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The matrix product structure is considered on a regular lattice in the hyperbolic plane. The 
phase transition of the Ising model is observed on the hyperbolic (5, 4) lattice by means of the 
corner-transfer-matrix renormalization group (CTMRG) method. Calculated correlation length 
is always finite even at the transition temperature, where mean-field like behavior is observed. 
The entanglement entropy is also always finite. 
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1. Introduction 

Classification of phase transitions is one of the central 
issue in the study of lattice models in statistical mechan- 
ics. When a system exhibits the second-order transition, 
normally the correlation length diverges at the transition 
point. As a result of scale invariance at the criticality, the 
transitions are characterized by scaling indices, where 
their values are completely classified in two-dimension 
by means of the conformal field theory. 

The mean-field like 2nd-order transition is exceptional 
in the point that the correlation length may not play an 
important role, in particular when the transition is de- 
scribed by the Landau free energy that is expressed as a 
simple polynomial of the order parameter.^' In this ar- 
ticle we focus on the mean-field like transition observed 
for the Ising model on the hyperbolic lattices, ^"^^^ the 
regular lattice in two-dimensional (2D) plane with con- 
stant negative curvature. ^^^ Among the hyperbolic (p, q)- 
lattices, which are the tessellations of the regular p-gons 
with the coordination number q,^'^' we consider the (5, 4)- 
lattice shown in Fig. 1 as an example. We calculate the 
correlation length ^ and entanglement entropy S in the 
neighborhood of the second-order transition temperature 
Tg , and judge whether or not the system is critical at this 
temperature. 

In the next section we explain the matrix product 
structure of the Ising model on the (5, 4)-lattice. We em- 
ploy the corner transfer matrix renormalization group 
(CTMRG) method, i-*'!^) a variant of the density ma- 



trix renormalization group (DMRG) method^^"^^' ap- 
plied to 2D classical models, ^^-' to obtain the thermody- 
namic properties of the model. We show the calculated 
results on ^ and S in §3. Conclusions are summarized in 
the last section. 

2. Matrix Product Structure on the Hyperbolic 
Lattice 

Consider the ferromagnetic Ising model on the (5,4)- 
lattice shown in Fig. 1. Each pair of neighboring sites is 
on a geodesic, which is drawn either by a line that passes 
through the center of the disk or by an arc. When there 




Fig. 1. The hyperbolic (5,4) lattice drawn in the Poincare disk. 
Open circles denote lattice points, where there are Ising spin 
variables. Those regions denoted by P correspond to half-column 
transfer matrices. 



is no external magnetic field, the Hamiltonian is given 

by 



H 






(2.1) 



where CTj = ±1 denotes the Ising spin variable at the z-th 
site, and where J > is the coupling strength between 
neighboring pair of sites denoted by {ij). It is convenient 
to introduce the interaction-round- a- face (IRF) Boltz- 
mann weight 

W,,Mm = (2.2) 



exp 






'k'^i 



for each pentagon, where i, j, fc, i and m denote the sites 
around it, and where /3 represents the inverse temper- 
ature. The partition function of the system is then ex- 
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pressed as 



config. (ijkim) 



(2.3) 



where the product is taken over all the pentagons de- 
noted by the group of sites {ijk£m), and where the sum 
is taken over all the possible spin configurations. We as- 
sume that the system is sufficiently large, and we con- 
sider the thermal equilibrium state deep inside the sys- 
tem, where the part is far from the system boundary. 

Let us observe the structure of the (5,4)-lattice. The 
thick straight (horizontal) line drawn in Fig. 1 divides 
the whole lattice into upper and lower halves. Let us 
introduce the new labeling 



^(0) (0) (0) (0) 



(2.4) 



from left to right for those spins on this line. We use 
the notation a^^' when we refer to the whole part of this 
horizontal row of spins. 

We next observe the thick arc that is perpendicular to 
the horizontal line we have considered, and that passes 
through the ^-th site a\ . We label those spins on this 
vertical geodesic as 



(-1) (0) (1) (2) 
, u fi . Of, , Of, . O I) . 



(2.5) 



from downward to upward. We use the notation erf when 
we refer to the whole part of this vertical column of spins. 
For the latter convenience, we introduce the half-infinite 
spin columns 



(-3) (-2) (-1) 
. . . , (7^ ,0-^ ,0-^ 

(1) (2) (3) 
Of, , o /) . u I) .... 



(2.6) 



for the lower and the upper part of the column spin ai' 
It should be noted that the spins 



.,0- 



0) Jo) Jo) Jo) 



e+i''^e ' 



-1' "£-2' ■ 



(2.7) 



for j =^ arc not on a geodesic, and therefore they cannot 
be regarded as a row spin. We have not explicitly shown 
the system boundary, and have not given any special la- 
bel to the boundary spins. When the system is off critical, 
we do not have to consider about the system boundary 
so strictly. 

We introduce the column-to-column transfer matrix 



^£+l,<? = ^('^£+lk£"): 



(2.8) 



which is created by multiplying all the IRF weights inside 
the stripes between £ + 1-th and l-th vertical geodesies, 
and taking spin configuration sum except for those col- 
umn spins o'lv^ and af. According to the division of the 
column spin 



^D JO) ^U 



(2.9) 



we can also express the column to column transfer matrix 
as a product of its upper and lower parts 

^ — T/'^D J^) ^U |„D ^(0) U^ 

= Pi+ijPi+i,i-> (2-10) 

where P^i £ and P^i £ are half-column transfer matrices 



(HCTMs) 

pD _ PC D (0) I D JO). 

Pe\e - P(af^\,<i|af,^U)^ (2.11) 

Using these notations the lower and the upper halves 
of the system are represented by the IRF-type matrix 
products^") 

config. i 
D 

config. i 
U 

(2.12) 

respectively, where configuration sums are taken for all 
the aP and a^ . Note that the partition function Z in 
Eq. (2.3) is expressed as the inner product between ^^ 
and *u 



Z = ^VI/D(a("))vl/U(^(0))^ 



(2.13) 



where the configuration sum is taken over the horizontal 
row spins. 

Statistical property of the left half of the system can 
be represented by the successive product of the column- 
to-column transfer matrices 



$^a„D,ar,a„U)=vL [] ^^+i. 



(2.14) 



,e>o 



where V'^ is a vector that specifies the boundary condi- 
tion at the left border of the system. In the same manner 
the property of the right half of the system is represented 
as 



$Vo^4"\-o^)= n^«-iW' 



(2.15) 



ve<o 



We assume that a very weak symmetry breaking field is 
imposed to the boundary spins, though we do not ex- 
plicitly refer to this condition in the following. Since the 
lattice structure shown in Fig. 1 is isotropic, it is also 
possible to represent $'^(ct(j') and ^^{ctq) as products of 
the half-row transfer matrices in the same manner as we 
have expressed *°(cr(°)) and *"(cr(°)) in Eq. (2.12).^!) 
Taking a partial contraction between $'"(cr5') and 
$'^(crQ'), we obtain the density matrix 

(2.16) 



D /-Di D\ 

P = P[<^0 FO ) 



Y: $^ao^ar,ao^)$^(ao^ar,a^ 



that is treated in the context of the bipartite entangle- 
ment to the vertical direction, or the block diagonal den- 
sity matrix 



/-D (0)| D (0)^ 



(2.17) 
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Fig. 2. Correlation Length ^ with respect to the temperature T. 



Fig. 3. Entanglement entropy of the MPS in Eq. (2.12). 



that is the 4-th power of the corner transfer matrix 

(CTM).l' 22-24) 

3. Correlation Length Obtained by CTMRG 

The matrix product representation of '^^{a^^') and 
^u^iy(o)^ in Eq. (2.12) has the same form as those 
for the square lattice Ising modeL Thus we can apply 
Baxter's variational formulation^' 22-24) qj. ^Yiq DMRG 
method^^"^^^ for the calculation of the free energy and 
other thermodynamic functions. In this article we em- 
ploy the CTMRG method, i". i5) a variant of the DMRG 
method, to obtain the correlation length ^ and the en- 
tangle entropy 5* of the Ising model on the (5, 4) lattice. 

The CTMRG method maps the half-column spins of 
and o-y , respectively, to block spins Qf ^^^^ CY ^y means 
of the the renormalization group (RG) transformation, 
which is obtained by the diagonalization of the density 
matrix in Eq. (2.17), or that in Eq. (2.16). Through 
this RG transformation, the HCTMs are mapped to the 
renormalized ones 



P^ 



D 

1+1,1 



— ^\^e+iT"e+i\^e ^"e ) 



M; 



P> 



u 



.(") /'UN 



£+1,1 — -P('''£+l'C£+ll "■£ jQ)j (3-1) 



where we put """ marks on top of renormalized matrices. 
The column-to-column transfer matrix is renormalized in 
the same manner 

f(y-B (0) >U iz-D (0) .UN 



T, 



e+i,e 



P, 



D 

1+1,1 -^e+i,e 



P, 



u 



(3.2) 



In the previous studies we have shown that the system 
exhibits the mean-field like second-order phase transi- 
tion, ^^^^^ where the spontaneous magnetization 



E(0) />D (0)1 aD (0)n 



M 



,D^(0) 



to '^0 



(3.3) 



below the transition temperature Tg is proportional to 



Figure 2 shows the correlation length 



e = 



1 



log Aq - log a 1 



(3.4) 



calculated from the largest eigenvalue Aq and the second 
largest one A^ of 2^,^,^^. We have regarded the interaction 
parameter J as the unit of energy, and set the lattice 
constant as the unit of length. We keep at most to = 40 
states for block spins, and actually to, = 5 is sufficient 
enough to draw the figure. It is clear that ^ is of the 
order of the lattice constant even at the second order 
transition point Tq = 2.799. The fact shows that the 
system is always off-critical, and thus the transition point 
cannot be called as the critical point. The entanglement 
entropy 



S = -E^ilog'^j 



(3.5) 



where u^ 
matrix fP 



is the z-th eigenvalue of the reduced density 
= P^iColCa) in Eq. (2.16), shown in Fig. 3 
is also finite for any temperature T. The fact coincides 
that a very small number of states to- = 5 is sufficient for 
getting thermodynamic quantities precisely. ^^^^-^ 

4. Conclusions and Discussions 

We have calculated the correlation function ^ and the 
entanglement entropy S of the Ising model on the hy- 
perbolic (5,4) lattice. Both of them remains finite at the 
transition temperature Tq . Therefore the mean-field like 
second-order phase transition observed for this system is 
not related to critical phenomena with diverging ^. This 
off-critical behavior is common to phase transitions ob- 
served in the tensor product formulations, 2^-29) -vvj^ej-e 
the trial state is finitely correlated due to the state num- 
ber limitation for the block spin variables. 

The mean-field nature of the phase transition might 
be explained by the path integral representation of the 
Green function on the lattice, if dominant contribution 
comes from the shortest path between two points. It 
should be noted that in the hyperbolic plane any de- 
viation from the shortest path causes increase of path 
length more than that on the flat plane. Such an increase 
reduces entanglement between two points. It should be 
noted that the Hausdorff dimension of the (5,4) lattice, 
or more general the hyperbolic {p, q) lattices, is infinite. 

From the path integral picture on the hyperbolic plane. 
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one reaches a deformation to ID quantum Hamiltonians, 
so called the hyperbolic deformation. ■^"'^^^ We conjec- 
ture that mean-field behavior also appears in ground- 
state phase transitions of deformed ID quantum systems, 
when NX is sufficiently larger than unity, where N is the 
system size. 

It is known that Kasteleyn formalism is applicable for 
dimer models on the hyperbolic lattices. ■^^•' Thus there 
is a mathematical interest on the hyperbolic lattices to 
find out a commutable row-to-row or column-to-column 
transfer matrices. 
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